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摘 要 :在 非 线 性 结构 的 振动 控制 设计 中 结构 模型 的 阶 数 不 宜 过 高 ,为 此 ,本 研究 以 含 局 部 非 线 性 

的 悬臂 梁 为 研究 对 象 , 开 展 影响 POD 降 阶 方法 所 得 低 阶 模型 精度 的 研究 。 着 重 分 析 了 非 线 性 强 

弱 、 降 阶 模 型 的 阶 数 、POD 模 态 获取 源 信号 的 激 振 类 型 响应 信号 的 采样 频率 和 响应 信号 采样 时 长 
等 因素 对 降 阶 模型 响应 预测 精度 的 影响 。 结 果 表 明 :对 于 强 非 线 性 的 局 部 非 线性 悬 辟 梁 系统 ,POD 

方法 同样 适用 ;在 选取 源 信 号 的 激 振 类 型 时 ,应 避免 选取 脉冲 激励 信号 ;响应 的 采样 频率 与 时 长 不 

一 定 要 选取 过 大 。 最 后 ,提出 了 一 种 针对 含有 噪声 信号 应 用 POD 方法 的 解决 方案 ,可 为 工程 应 用 
= 提供 有 益 的 参考 。 
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Adaptability analysis of reduction model of a nonlinear 
beam based on POD method 
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Abstract :In the vibration control design of nonlinear structures, the order of the structure model should 
not be too high. For this reason , this paper takes the cantilever beam with local nonlinearity as the research 
object to carry out the research on the low-order model accuracy obtained by the POD reduction method. 
The analysis focuses on the influence of factors such as the strength of nonlinearity , the order of the re- 
duced-order model,the excitation type of the source signal obtained by the POD mode, the sampling fre- 
quency of the response signal and the sampling duration of the response signal and other factors on the re- 
sponse prediction accuracy of the reduced-order model. For strongly nonlinear local nonlinear cantilever 
beam systems, the POD method is also applicable. When selecting the excitation type of the source signal, 
it is necessary to avoid the conclusion that the sampling frequency and duration of the pulse excitation sig- 
nal and response selected should not be too large. Finally ,a solution for applying the POD method to noisy 


signals is proposed , which can provide a useful reference for engineering applications. 
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近年 来 ,航天 领域 由 轻薄 材料 制 成 的 天 线 和 大 
能 板 的 振动 控制 问题 是 人 们 的 研究 热点 ;在 机 
a a t 


土木 工程 领域 由 于 风 载 和 地 震 引 起 的 大 型 结构 振动 
问题 也 日 益 突 出 ”1 。 过 去 的 几 十 年 里 有 限 元 理论 、 


方法 及 其 应 用 取得 了 巨大 的 进步 ,对 于 静 力 分 析 , 大 
规模 计算 模型 的 计算 效率 还 可 以 接受 ,但 结构 动力 
分 析 的 计算 量 要 比 静 力 分 析 大 几 个 量 阶 , 几 乎 不 可 
能 采用 这 样 的 模型 进行 结构 动力 分 析 , 并 且 由 于 控 
制 系统 设计 的 复杂 性 ,控制 模型 的 规模 不 能 太 高 ,如 
果 结 构 系统 的 自由 度 太 大 ,就 会 产生 控制 器 及 观测 
器 < 洲 出 ”等 问题 。 

全 许多 结构 并 非 由 单独 一 根 梁 单元 组 成 ,而 是 由 
2 梁 单 元 组 成 ,所 以 对 连接 结构 的 研究 也 是 人 们 
注 的 一 大 问题 。 目 前 也 有 大 量 的 研究 通过 非 线 
性 改作 为 系统 基础 线性 模型 中 的 内 部 反馈 力 ,并 通 
过 实验 进行 识别 非 线性 系数 , 孙 伟 等 "建立 了 
螺栓 连接 梁 的 非 线 性 解析 模型 ,并 通过 实验 辨别 相 
关 允 数 , 在 弯曲 方向 上 可 将 螺栓 连接 结构 看 成 是 包 
食 记 方刚 度 项 的 形式 ,所 以 本 研究 以 含 局 部 非 线性 
的 沪 避 梁 为 模型 进行 研究 。 

仿 本 征 正 交 分 解 ( proper orthogonal decomposition , 
pb 吕方 法 是 一 种 常见 的 基于 投影 In EET E, 
HEE POD 方法 得 到 的 降 阶 模型 自由 度 较 低 , 且 为 原 
系统 的 最 佳 平方 通 近 在 高 维 复杂 系统 计算 效率 及 精 
度 和 方面 优势 显著 ,该 方法 也 已 广泛 用 于 许多 研究 
领域 ,以 减少 实际 工程 系统 的 阶 数 。 例 如 ,在 流体 力 
学 的 背景 下 ,KUNISCH 等 …” 分 析 了 将 POD 方法 用 
于 流体 力学 一 般 方程 。 在 结构 动力 学 领域 , KUMAR 
等 中 使 用 POD 方法 获得 非 线 性 结构 动力 系统 的 降 阶 
模型 。TABIB 等 ' ”使 用 POD 方法 来 设计 化 学 反应 器 
过 程 的 经 验 基础 函数 。 本 研究 选用 POD 方法 作为 基 
本 降 阶 方法 对 含 局 部 非 线 性 悬臂 梁 进 行 探 究 。 

本 研究 针对 含 局 部 非 线性 梁 的 模型 降 阶 问题 ， 
通过 误差 统计 的 方法 ,研究 了 应 用 POD 方法 时 , 非 
线性 强 弱 、 阶 数 、 激 振 方法 和 采样 方法 对 所 得 降 阶 模 
型 的 精度 影响 ,并 提出 一 种 基于 实际 工程 中 含 噪 声 
响应 信和 号 建立 降 阶 模型 的 方案 。 本 研究 得 到 了 含 局 
部 非 线 性 梁 的 低 阶 模型 ,对 进行 非 线性 结构 的 振动 
控制 设计 有 参考 意义 。 


1 局 部 非 线性 粱 的 有 限 元 模型 


本 研究 针对 含有 局 部 非 线性 特性 的 悬臂 梁 模 
型 ,如 图 1 所 示 , 其 局 部 非 线 性 由 梁 弹性 支承 端的 非 
线性 立方 弹簧 产生 。 线 性 梁 结 构 采 用 欧 拉 梁 模 型 建 
模 , 即 不 计 横 截面 转动 惯量 和 剪 切 变形 效应 。 


= 


图 1 含 局 部 非 线性 的 悬臂 梁 模型 示意 图 


Fig.1 Schematic of the nonlinear beam 


沿 轴 线 取 欧 拉 梁 微 段 进行 受 力 分 析 , 微 段 示意 
"aues 
pot) 
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图 2 m 


Fig.2 Micro segment of beam under stress 


通过 对 梁 微 段 受 力 分 析 后 ,可 得 到 欧 拉 梁 弯 曲 
振动 的 微分 方程 为 
了 
EI a -m(x) n = p(x,t) (1) 


其 中 : v(x,t) 2938 B 8 BE; 0(x,0) 为 截面 转角 ; 
p(x,t) 为 作用 在 梁 上 的 载荷 ; M(x,t) REUSI; 
V(x ,t) 为 截面 剪 力 。 
欧 拉 梁 弯 曲 振 动 单元 参数 如 图 3 所 示 。 
y 


图 3 欧 拉 梁 弯 曲 振动 单元 参数 


Fig.3 Parameter of Euler-beam flexural vibration 
根据 最 小 势能 原理 ,最 终 可 得 欧 拉 梁 有 限 元 单 
元 的 方程 
K'B*4g' M.'6.°=F. (2) 
K,* 为 欧 拉 梁 的 单元 刚度 矩阵 (对 称 矩 阵 ) , H. 
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6 31 -6 3l 
I 2P 350] P 
kaea s ? (3) 
/|-6 -31 6 9M 
31 P -3 27 


M“ 为 欧 拉 梁 的 单元 质量 矩阵 ( 对 称 和 矩阵 ) , 且 


156 221 54 -13l 
2 2 
M. - Apl| 221 4l 131 3l (4) 
420| 54 13] 156 -22l 
-131 -3P -21 4P 


Hop: /为 梁 长 ;上 为 梁 的 杨 氏 模 量 ; I WERE; 4 
为 梁 界 面 面 积 ; p 为 梁 的 密度 。 

根据 位 移 协调 条 件 和 节点 平衡 条 件 , 可 组 装 得 
到 梁 总 体 刚度 矩阵 和 总 体质 量 抢 阵 ,并 加 入 悬臂 梁 
的 边界 条 件 ,可 得 到 欧 拉 染 的 有 限 元 动力 学 方程 , 即 


: Mž+Kx=0 (5) 
孔 假 设 该 非 线性 梁 中 的 阻尼 形式 为 瑞 利 阻 尼 
O C = aM + BK (6) 
HISP TES ep SERO RU vr 77 9896 A8 JE AR PEOR 
N(x) =- kx (1) 


JR 为 非 线 性 系数 。 最 终 得 到 离散 非 线 性 梁 的 
373^ 71 FON 


= Mž + Cx + Kx + N(x) =F (8) 
H, F 为 施加 于 梁 上 的 外 激励 ,其 具体 形式 依据 外 


激励 信号 确定 。 


3 F POD 方法 的 含 局 部 非 线 性 梁 
三 的 降 阶 模型 


“对 于 高 维 复杂 系统 (n 维 ) 的 位 形 空间 变量 y(x， 
t) ,相应 数值 仿真 信号 或 实验 信号 的 离散 时 间 序 列 
可 表示 为 y(x,t;) ,i=1,…,N, 可 将 其 投影 到 一 组 完 
备 正 交 规范 基 o = (e, (x) 及 -1 张 成 的 空间 上 , 即 

y(x,t,) = 2 aye GO) i -]1,-—,N (9) 
其 中 和 表示 某 一 时 刻 的 数据 索引 ; wV 表示 离散 数据 
总 数 尖 表示 某 一 个 正 交 规范 基 索 引 ; 表示 正 交 规 
范 基 总 数 , 也 为 空间 维 数 。 

寻找 出 该 正 交 规范 基 过 程 可 表示 为 满足 下 式 的 
最 小 约束 问题 。 

min( [y(x,5) -Py °), lel? =1 (10) 


KAREA INVOER 425 


该 问题 可 通过 Lagrange 乘 子 法 求解 。 
Jle] = <10) I -Alel2-D (1) 
极 值 的 必要 条 件 为 泛 函 为 0, 即 5J1 =0。 因 此 ， 
通过 求解 以 下 公式 的 本 征 函 数 ,可 以 得 到 正 交 基 


Í, (y(x)y(x))e(x)dx = Ap(r) (12) 


上 式 为 第 二 类 Fredholm 7j fé, HAX RAON A 4H 
关 函 数 : RCx,x) =《〈y(x)y(x )〉, 有 限 维 情况 下 ， 
POD 降 阶 模 态 为 下 属 自 相关 和 矩阵 R 的 特征 向 量 。 


1 
R = — [y], (13) 


N, 
其 中 ,NN, 为 采样 点 数 。 
将 高 维 位 形 空间 变量 投影 到 前 少数 几 个 POD 
模 态 张 成 的 子 空间 上 (特征 向 量 按 对 应 特征 值 降 序 
排列 ) , 则 能 以 最 少 的 降 阶 模 态 数 m, 获 得 高 维 复杂 
系统 位 形 空间 变量 在 采样 时 间 内 的 最 佳 平方 逼近 。 


y(x,t;) 一 Y uu Ge) i -2]1,--,N;m-«n 
fi 


(14) 

定义 坐标 y = Pu, P 为 正 交 规范 基 gp 张 成 的 子 

ZH), P = [P1,925 Pn]nxm 4 为 POD 模 态 坐标 ， 

u(t) = [u (t) us) sus (1) ]" o 将 其 带 入 非 线 

性 梁 的 有 限 元 模型 的 动力 学 方程 式 (8) 中 ,化 简 
可 得 

Mii+Kiu+Cu+N(u) =F (15) 

其 中 : M = P'MP,C = P'CP,K =P'KP,N(u) = 


PTIN(x) ,F =P'F, 分 别 为 降 阶 模型 的 质量 .阻尼 、 刚 
度 矩 阵 及 非 线 性 力 向 量 和 外 载荷 力 向 量 。 由 此 可 得 
含 局 部 非 线 性 梁 的 降 阶 模型 。 


3 含 局 部 非 线性 梁 降 阶 模型 响应 精度 
分 析 


将 含 局 部 非 线 性 的 梁 离散 为 10 个 单元 , 则 xs 
RR”, 从 而 得 到 40 维 的 相 空 间 。 依 据 参 考 文 献 [20] 
中 梁 模 型 , 梁 的 几何 和 材料 参数 如 表 1 所 示 。 
通过 求解 该 梁 结构 的 特征 方程 组 
IK-wM|I=0 (16) 


其 中 :〈:…〉 表 示 采 样 时 间 段 内 的 平均 算 子 ;P 表示 投 
影 算 子 ; Cu 表示 Hilbert 空间 上 的 内 积 ; ‖ .| 
den LD 的 范 数 。 


可 以 得 到 该 梁 的 固有 频率 ,前 4 阶 结 果 如 表 2 所 示 。 
本 研究 所 有 计算 采用 的 是 自 编 有 限 元 程序 ,为 
证 明 结果 可 靠 性 ,在 相同 参数 和 初始 条 件 下 与 商业 
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软件 ABAQUS 计算 结果 对 比如 图 4 所 示 。 
R1 ， 梁 的 几何 和 材料 参数 


Tab.1 Geometric and material parameters for the beam 
符号 值 单位 
l 2 700 mm 
A 100 mm? 
p 1780 x107? kg/mm 
E 45 x 10° kPa 
I 833.3 mm* 
K 6 N/mm? 
a 1.25 x10 * s! 
B 2.5 x10^* s 


R2 梁 的 固有 频率 
Tab.2 Natural frequency for the beam 


阶 数 1 2 3 4 
推导 值 /Hz 7.0005 43.8728 122.8724 240.9489 
"ABAQUS 值 /Hz | 6.9994 43.8500 122.7600 240.5400 
Toc 数值 结果 
一 一 ABAQUS 仿 真 


:/(10* m) 


x. 


Tg 


-4 
0 50 100 
时 间 /s 


图 4 自 编程 序 与 ABAQUS 结果 对 比 


Fig.4 Comparison with ABAQUS results 
从 图 4 中 可 以 看 出 自 编程 序数 值 仿真 与 
ABAQUS 仿真 结果 相同 ,通过 模 态 分 析 , 可 得 到 与 自 
编程 序 相同 的 固有 频率 ,证 明了 自 编程 序 的 可 靠 性 。 
曲线 性 梁 的 频 域 响应 如 图 5 所 示 , 取 [1 Hz, 
50 Hz] 频率 范围 ,覆盖 前 两 阶 固有 振动 频率 。 


幅 

| 
[eN 
© 


-100 
0 10 20 30 40 50 
频率 /Hz 
图 5 线性 梁 的 频 域 响 应 
Fig.5 The frequency response of the linear beam 
在 梁 的 端 部 施加 不 同 幅 值 的 激 振 力 , 梁 在 一 阶 
固有 频率 附近 的 响应 呈现 非 线 性 特性 。 


ee 1.05 x 10? N 非 线性 


A 一 -一 1.00X10 NIERE 
E = - 095X103N 非 线性 


-10 —— 1.00 X 10° N&R TE 


ea 
id 
gj -15 
EE 

-20 

E25 

-30 

6.90 6.95 7.00 7.05 7.10 
频率 /Hz 


图 6 梁 的 非 线性 特性 
Fig.6 The nonlinear characteristics of the beam 

图 6 中 实 线 为 线性 梁 幅 频 特性 ,虚线 为 施加 非 
线性 弹簧 后 的 幅 频 特性 ,从 图 中 可 以 看 出 施加 非 线 
性 弹簧 后 会 表现 出 的 刚度 硬化 现象 。 

为 研究 POD 降 阶 方法 对 含 局 部 非 线 性 粱 的 适 
应 性 ,引入 误差 值 判断 降 阶 模型 的 响应 精度 ,定义 

平均 绝对 误差 


1 n 
A 2 |y; — Yiron) | (17) 
平均 相对 误差 
8 = LY |n -roo (18) 
n izi Yi 


其 中 : 为 求解 过 程 中 的 时 间 节 点 总 数 ; y 为 在 时 间 
点 i 处 的 未 降 阶 模型 响应 结果 ; yop, 为 降 阶 模型 在 
时 间 点 i 处 的 响应 结果 。 

本 研究 中 所 有 误差 均 为 谐 波 激励 下 原 模 型 响应 
与 POD 降 阶 模型 响应 的 误差 , 且 在 该 模型 算 例 中 可 
接受 的 误差 最 大 值 为 10% 。 


3.1 非 线性 强 弱 对 降 阶 模型 预测 精度 的 影响 


POD 降 阶 方法 主要 是 在 线性 系统 的 基础 之 上 提 
出 的 一 种 降 阶 方法 ,所 以 对 其 在 非 线性 系统 中 的 适 
应 性 问题 的 探究 是 本 研究 的 研究 基础 。 下 面 讨 论 含 
非 线 性 梁 模 型 中 不 同 非 线性 系数 c 对 预测 精度 的 
影响 。 

在 非 线性 悬臂 梁 的 弹性 支承 端 部 施加 谐 波 激励 
F = F,cos(ot) , Fo BU1 x 107N , 激 振 频率 ww = 
30 Hz ,为 排除 由 于 采样 频率 过 低 和 采样 时 长 过 短 所 
产生 的 误差 精度 影响 ,采样 频率 取 1 000 Hz ,采样 时 
长 取 1000s , 则 在 该 参数 下 可 得 到 系统 的 响应 信号 
y = [nasa Y2 l2x100000 。 以 该 信号 为 源 信 
号 ,构造 POD 降 阶 模 态 。 以 梁 弹 性 支承 端点 为 分 析 
点 ,对 比 未 降 阶 系统 和 通过 POD 降 阶 的 系统 在 上 述 
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参数 激励 下 的 结果 。 

在 梁 的 端 部 施加 相同 幅 值 的 激 振 力 和 不 同 的 非 
线性 大 小 , 则 梁 在 一 阶 共振 激励 频率 附近 的 幅 频 响 
应 特性 如 图 7 所 示 。 
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图 7 不 同 非 线 性 系数 下 的 幅 频 响应 

Fig.7 The nonlinear characteristics 
三 取 3 阶 进 行 POD 降 阶 模 态 构造 , 非 线性 系数 K 

ZAH 6 .600 .6 000 ,针对 这 3 种 非 线性 系数 的 未 降 
阶 天 系统 模拟 数据 ,分 别提 取 构 造 了 3 个 POD 降 阶 
模型 ,然后 分 别 与 原 系统 对 比 ,将 降 阶 后 模型 与 原 未 
降 阶 模型 分 别 在 激励 频率 为 1 ~ 50 Hz 的 响应 误差 
AR 比如 图 8 所 示 。 
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FT 图 8 非 线性 强 弱 误差 对 比 


Fig.8 Nonlinearity error comparison 
从 图 8 可 以 看 出 , 当 非 线性 系数 k 分 别 取 6、 
600 .6 000 时 , 仍 能 够 保证 平均 相对 误差 频 域 上 的 平 
均 误 差分 别 为 5% 1% 、3% 的 降 阶 精度 ,所 以 POD 
降 阶 方法 可 应 用 于 不 同 非 线性 强 弱 系统 。 
后 文 分 析 所 采用 的 几何 和 材料 参数 均 如 表 1 所 
示 , 即 非 线性 刚度 系数 为 6。 


3.2 模型 阶 数 对 预测 精度 的 影响 


POD 方法 是 在 保证 一 定 精 度 的 情况 下 ,将 高 阶 
系统 降 阶 为 低 阶 系统 的 一 种 方法 ,进而 能 够 大 大 节 
约 结构 响应 的 计算 时 间 。 下 面 讨 论 模 型 阶 数 对 预测 
精度 的 影响 。 

在 非 线性 悬臂 粱 的 弹性 支承 端 部 施加 与 上 节 相 
同 的 谐 波 激励 ,并 取 相 同 的 采样 方法 ,可 得 响应 信 
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号 。 以 该 信号 为 源 信号 ,构造 POD 降 阶 模 态 。 以 梁 
弹性 支承 端点 为 分 析 点 ,对 比 未 降 阶 系统 和 通过 
POD 降 阶 的 系统 在 上 述 参数 激励 下 的 结果 ,误差 对 
比如 图 9 所 示 。 
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*|9 阶 数 影响 
Fig.9 The order influence 

从 图 9 中 可 以 明显 看 出 ,在 全 阶 为 20 阶 的 系统 
中 , 取 3 阶 POD 降 阶 模型 即 可 得 到 0. 9096 的 相对 误 
差 。 本 研究 所 有 计算 均 为 自 编 程序 ,采用 MATALB 
软件 中 ode45 函数 进行 求解 。 表 3 列举 了 未 降 阶 模 
型 与 3 阶 POD 降 阶 模型 的 计算 时 间 。 可 以 看 出 , 降 
阶 模型 的 计算 时 间 仅 为 未 降 阶 模型 的 3% ,所 以 基 
于 POD 方法 所 得 降 阶 模型 是 较为 高 效 的 。 

表 3 计算 时 间 对 比 


Tab.3 | Comparison of calculation time 


项 目 未 降 阶 
计算 时 间 /s 


3 阶 POD 


281.56 8.50 


3.3 ”基于 不 同 激 振 响 应 信号 的 降 阶 模型 精度 分 析 
基于 不 同 谐 波 频率 激励 信号 的 降 阶 模型 


考虑 对 包含 前 两 阶 固 有 频率 的 预测 精度 的 影 
响 , 故 施加 谐 波 激励 时 , 激 振 频 率 也 应 包含 前 两 阶 
有 频率 。 在 非 线性 悬臂 梁 的 弹性 支承 端 部 施加 3.1 
节 相 同 的 谐 波 激励 , 激 振 频率 从 1 Hz 开始 ,以 1 Hz 
间隔 取 到 50 Hz ,采样 频率 取 1 000 Hz ,采样 时 长 取 
1000 su 从 3.2 节 可 知 , 当 选取 3 阶 进行 POD 降 阶 模 
态 构 造 时 , 便 可 达到 较为 合适 的 预测 精度 ,所 以 阶 数 
选取 3 阶 。 以 梁 弹 性 支承 端点 为 分 析 点 , 降 阶 结 
如 图 10 所 示 。 

在 图 10(b) 和 (ec) 中 ,采样 数据 频率 表示 获得 降 
阶 模型 时 ,施加 谐 波 激 励 响应 的 频率 ;实验 数据 频率 


表示 对 得 到 的 降 阶 模型 施加 谐 波 激 振 的 频率 。 
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(d) 30 Hz 处 位 移 时 程 图 (以 30 Hz 为 源 数 据 ) 
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Fig. 10 The excitation frequency influence 


从 图 10 误差 图 中 可 以 看 出 ,误差 精度 较 差 部 分 
存在 于 共振 峰 或 者 反共 振 峰 附近 , 即 在 该 梁 的 7、 
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31,44 Hz 附近 。 并 且 图 中 也 明显 可 以 看 出 , 若 要 得 
到 某 个 降 阶 响应 最 好 的 效果 ,应 从 该 频率 附近 的 激 
振 响 应 中 获取 POD 降 阶 模型 ,如 ,从 30 Hz 的 激 振 响 
应 中 获取 的 POD 降 阶 模型 在 施加 30 Hz 谐 波 激 振 频 
率 的 降 阶 模型 得 到 最 小 误差 。 但 从 整体 结果 看 , 均 
可 达到 相对 误差 小 于 5% 的 水 平 ,说 明 谐 波 激 励 响 
应 为 一 种 可 选 的 降 阶 信号 。 
3.3.2 基于 扫 频 激励 响应 信号 的 降 阶 模型 精度 

当 响 应 数据 中 含有 较为 丰富 的 频率 信息 时 ,可 
能 会 获得 较 好 的 预测 精度 , 而 扫 频 激励 信号 中 含有 
丰富 的 频率 信息 。 为 进行 验证 ,在 非 线性 悬臂 梁 的 
弹性 支承 端 部 采用 线性 扫 频 激励 , 即 w = wo vt, Dt] 


最 终 扫 频 信号 为 = Fosin(wot + Lon) ,其 中 F, 


取 1 x 107N , w WO Hz, v 取 0.05, 采 样 频 率 取 
1 000 Hz ,采样 时 长 取 1000s , 阶 数 选取 3 阶 。 以 梁 


弹性 支承 端点 为 分 析 点 , 扫 频 降 阶 结果 如 图 11 
所 示 。 
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Fig. 11 The sweep excitation influence 


MELL 可 以 看 出 ,从 扫 频 激励 响应 获得 的 降 阶 
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模型 的 误差 最 高 仅 为 0.21% ,说 明 扫 频 激励 响应 也 
为 一 种 可 选 的 降 阶 信和 号。 
3.3.3 脉冲 激励 对 预测 精度 的 影响 

由 于 阻尼 的 存在 , 当 对 该 非 线 性 梁 受 到 脉冲 激 
励 , 如 锤 击 时 ,其 响应 振幅 与 振动 频率 将 快速 衰减 ， 
所 以 尽管 脉冲 激励 的 响应 也 包含 着 丰富 的 频率 特 
性 ,但 由 于 高 阶 频率 衰减 过 快 ,会 对 降 阶 模型 精度 产 
生 不 利 影 响 。 在 非 线性 悬臂 梁 的 弹性 支承 端 部 施加 
脉冲 信号 = Posin((Cmi)《r),7 取 0.05 ,在 起 始 状 
态 下 (t<0.05 s) 施 加 载荷 Fe =0. 1 N ,之 后 染 将 处 于 
自由 振动 状态 ,采样 频率 取 1 000 Hz ,采样 时 长 取 
1000 s, 阶 数 选 取 3 阶 。 以 梁 弹 性 支承 端点 为 分 析 
点 ,结果 如 图 12 所 示 。 
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图 12 脉冲 激励 影响 


Fig.12 The pulse excitation influence 
从 图 12 可 以 看 出 ,从 脉冲 激励 响应 获得 的 降 阶 
模型 ,误差 最 高 能 够 达到 202% , 达 不 到 精度 要 求 ， 
所 以 在 应 用 降 阶 响应 预测 时 ,应 避免 从 基于 脉冲 激 
励 响 应 的 信号 中 获取 降 阶 模型 。 
3.3.4 基于 随机 激励 响应 信号 的 降 阶 模型 精度 
在 非 线性 悬臂 梁 的 支承 端 部 施加 随机 激励 , 假 
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UE 


设 随 机 激励 力 服从 高 斯 白 噪声 分 布 , 则 可 采用 平稳 
高 斯 过 程 的 三 角 级 数 法 近似 模拟 随机 激励 , 即 激 


WRR U) = Y, asinado) 。 取 噪声 振动 


强度 为 1 x10N, 限 带 白 噪声 频率 范围 为 0 ~ 
50 Hz ,频率 分 辩 率 为 0.25 Hz ,可 确定 模拟 函数 的 所 
有 参数 。 采 样 频率 取 1 000 Hz ,采样 时 长 取 1000 s, 
阶 数 选取 3 阶 。 以 梁 弹 性 支承 端点 为 分 析 点 , 结 

如 图 13 所 示 。 
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(c) 30 Hz 处 位 移 时 程 图 


图 13 ”随机 激励 影响 
Fig.13 The random excitation influence 

从 图 13 可 以 看 出 ,从 随机 激励 响应 获得 的 降 阶 
模型 ,误差 最 高 也 仅 为 0.21% ,所 以 可 从 扫 频 激励 
响应 中 获取 降 阶 模型 。 

在 相同 的 激励 幅 值 ,相同 的 采样 时 长 和 相同 的 
采样 频率 下 ,由 不 同 激励 响应 得 到 的 隆 阶 响应 预测 
对 比 结果 如 图 14 所 示 。 

结果 表明 ,基于 POD 方法 , 扫 频 激励 ,随机 激励 
和 谐 波 激励 较 脉冲 激励 的 响应 均 能 得 到 符合 精度 的 
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降 阶 模型 ,并且 特 别 在 使 用 激励 响应 时 ,应 优先 使 用 
降 阶 模型 所 应 用 频率 段 的 谐 波 激励 响应 。 
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3.4. 对 降 阶 模型 选取 合适 的 采样 频率 


X 


(0 高 频 采 样 频 率 所 获得 的 数据 信号 能 够 包含 高 频 
的 数据 频率 信息 ,进而 可 获得 更 好 预测 精度 ,这 是 人 
们 括 一 般 认 知 的 。 若 采用 低频 的 采样 频率 ,如 要 分 
HY 50 Hz 的 频率 信号 ,采用 100 Hz 的 采样 频率 ,是否 
还 能 得 到 较 好 的 预测 精度 。 对 于 相同 的 采样 时 长 ， 
低频 的 采样 频率 能 够 得 到 更 少 的 数据 长 度 ,更 有 利 
于 对 数据 进行 分 析 处 理 。 

在 非 线性 悬臂 梁 弹 性 支承 端 部 施加 与 3.3.4 节 
相同 的 随机 激励 ,采样 频率 取 100 Hz ,采样 时 长 取 
1 000 s, 阶 数 选取 3 阶 。 以 梁 弹 性 支承 端点 为 分 析 
点 ,与 采样 频率 为 1 000 Hz 的 对 比 结果 如 图 15 
所 示 。 

从 图 15 可 以 看 出 ,将 采样 频率 降低 到 极限 最 低 
分 辨 频率 100 Hz, 降 阶 精度 仍然 与 1000 Hz 采样 频 
率 的 相同 ,所 以 在 进行 降 阶 采样 频率 选取 时 ,可 从 低 
至 高 进行 选取 来 构造 降 阶 模型 ,在 该 模型 下 可 选取 
100 Hz 的 采样 频率 。 
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图 15 采样 频率 影响 
Fig.15 The sampling frequency influence 


3.5 对 降 阶 模型 选取 合适 的 采样 时 长 


长 时 的 采样 时 间 能 够 使 得 采样 数据 中 包含 更 多 
的 重复 的 稳 态 信号 ,而 丰富 的 频率 信息 一 般 包含 于 
结构 的 瞬 态 过 程 。 所 以 选取 较 短 的 采样 时 长 ,不 仪 
能 够 减少 采样 数据 ,节约 计算 时 长 ,而 且 还 能 够 滤 去 
由 于 数据 过 多 而 产生 的 元 余 信息 。 

在 非 线性 悬臂 梁 弹 性 支承 端 部 施加 随机 激励 ， 
采样 频率 取 1 000 Hz, 采 样 时 长 取 50 s, 阶 数 选取 3 
阶 。 以 深 弹 性 支承 端点 为 分 析 点 ,与 采样 频率 为 
1 000 s 的 对 比 结果 如 图 16 所 示 。 
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(b) 平均 绝对 误差 对 比 (c) 平均 相对 误差 对 比 
图 16 采样 时 长 影响 
Fig.16 The sampling length influence 


从 图 16 中 可 以 得 出 ,采样 时 长 取 50s 与 1000s 
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所 得 到 的 降 阶 模型 的 误差 相差 非常 小 ,并 且 在 部 分 
频率 段 内 50 s 采样 时 长 得 到 的 降 阶 模型 的 误差 小 于 
1 000 s 采样 时 长 的 误差 ,说明 在 选取 采样 时 长 时 ,可 
选择 短 时 长 进行 降 阶 模型 ,在 该 模型 下 可 选取 50 s 
的 采样 时 长 。 


4 对 含 噪声 响应 信号 进行 降 阶 模型 分 析 


在 实际 的 振 劲 测量 中 ,由 于 传感器 的 选择 , 传 感 
髓 安装 方式 , 测 点 选取 和 被 测 对 象 的 工 况 环境 等 ,都 
会 使 得 所 测 得 的 振动 信号 中 伴随 大 量 噪声 ,所 以 研 
究 含 噪声 信号 对 降 阶 模型 的 精度 的 影响 在 实际 的 工 
程 应 用 中 有 重要 的 意义 。 
廊 对 响应 结果 施加 有 一 定 信 噪 比 的 噪声 , 信 噪 比 
ELUF. 
P, 


Z dB) = 10 log [ , 


ignal ) ES 20 logi ( Punt] 


noise noise 


OHER MER REE ARREDI 30 Hz 与 
3. 本 他 相 同 的 谐 波 激励 ,采样 频率 取 1000 Hz ,采样 
时 攻取 1000 s ,假设 响应 信号 中 含有 信 噪 比 为 10 dB 
的 响声 ,以 梁 弹 性 支承 端点 为 分 析 点 ,施加 噪声 的 信 


三 这 原 未 施加 噪声 信号 响应 结果 对 比如 图 17 所 示 。 
6 
dT. 一 全 me 
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频率 /Hz 
(b) 快速 传 里 叶 变 换 


图 17 不 同 响应 结果 对 比 
Fig.17 The response comparison 
基于 噪声 信号 进行 降 阶 分 析 , 结 果 表 明 , 加 入 噪 
声 后 对 预测 精度 有 非常 大 的 影响 。 考 虑 选择 合适 的 
降 噪 方法 进行 去 噪 后 进行 POD 降 阶 模 态 的 计算 。 


ade rcm pfe pe IV FRAT asi 


的 降 噪 方法 ,SVD 重 构 法 ,该 降 品 过 程 如 下 。 
1) f43& Toeplite ERE Y((N -元 +1) xL) 


xXL-1) y(L-2) y(0) 
ys X2 —»xL-D y(1) 
yN -1) y(N-2) yN - L) 


其 中 : y(n) 为 合 噪 声 的 信号 ,0 < n <N ; 工 为 预测 
阶 数 ,一 般 选择 二 = 3 。 
2 ) 计 算 含 噪 矩 阵 了 的 SVD 分 解 , 即 了 = LI 。 
3) 令 x, (n) 表示 使 用 秩 * 近似 合成 的 信号 ,并 
为 秩 ， 近似 定义 以 下 能 量 差 。 
E, Et (n)1- Ela} (n)1 
EXA) = |E, - oil, k-12,,L, RÆ 


REIR AC) , 则 可 估计 有 效 秩 r 。 
4) 计 算 了 的 秩 7 近似 = Y our! Eb u, 
RI v, 分 别 为 了 的 左右 奇异 向 量 。 


5) 对 成 的 对 角 线 元 素 进行 平均 运算 ,得 到 新 
矩阵 


x(L-1) x(L-2) x, (0) 
X, au dua 1) aeg 
(UN - 1) x (N -2) «(N - L) 


从 第 一 行 和 第 一 列 , 可 以 基于 秩 7 近似 的 合成 
信号 x,(n),n =0,1,2,…,N-1。 

本 研究 将 该 降 噪 方法 与 POD 降 阶 方法 结合 ,得 
出 含 噪声 响应 的 降 阶 模型 。 该 方法 的 关键 在 于 找到 
合适 的 有 效 秩 , 秩 选取 过 高 ,会 对 噪声 的 拟 合 更 加 精 
确 , 选 取 过 低 会 导致 无 法 还 原 原 数据 。 在 本 研究 中 ， 
该 有 效 秩 为 r=6。 

由 于 需要 构造 Toeplite 矩阵 ,所 以 较 大 的 数据 长 
度 在 执行 该 过 程 时 会 有 些 乏 力 ,但 由 3.4 节 与 3.5 
节 可 知 ,对 该 模型 ,可 选取 采样 频率 100 Hz ,采样 时 
K 50 s, 进 行 含 有 10 dB 的 噪声 的 谐 波 激 励 响应 信和 号 
进行 降 噪 ,对 比 结果 如 图 18 所 示 。 

从 图 18 可 以 看 出 ,该 方法 能 够 有 效 的 去 除 信和 号 
中 所 包含 的 噪声 。 以 降 噪 信号 为 源 信号 ,构造 POD 
降 阶 模 态 , 当 继 续 采样 阶 数 为 3 阶 时 ,发现 降 阶 模型 
精度 仍然 较 差 。 可 能 是 降 品 后 的 信号 所 包含 的 信息 
与 原 信号 仍 有 差距 ,但 继续 研究 发 现 , 提 高 阶 数 便 能 
够 提高 降 阶 精度 , 当 阶 数 选取 至 12 阶 时 , 可 得 到 符 


LOIZOUI2 介绍 了 一 种 对 含有 高 斯 白 噪 声 信 和 号 


合 标准 误差 的 结果 ,最 高 相对 误差 仅 为 2% 。 以 梁 
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弹性 支承 端点 为 分 析 点 ,对 比 结果 如 图 19 所 示 。 
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图 18 不 同 响应 结果 对 比 响应 


Fig. 18 The response comparison 
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19 ”误差 对 比 
Fig. 19 The error comparison 


本 节 对 含有 噪声 的 响应 信号 进行 讨论 , 发现 


SVD 重 构 法 可 以 高 效 的 滤 去 含有 高 斯 白 噪 声 的 信 
号 ,并 可 以 使 用 提高 阶 数 的 方法 提高 降 阶 精度 ,从 而 
获得 低 阶 的 POD 降 阶 模型 。 


论 


5 结 

本 研究 建立 了 含 局 部 非 线 性 梁 的 一 般 模型 ,从 
简单 模型 揭示 基于 POD 的 模型 降 阶 方法 的 适应 性 
问题 。 通 过 数值 分 析 ,得 到 以 下 几 点 结论 。 

1) 基 于 POD 的 模型 降 阶 方法 对 局 部 非 线性 梁 
动力 学 问题 具有 良好 适用 性 ,可 应 用 于 强 非 线性 系 
统 中 ,对 于 本 研究 的 研究 对 象 ,3 阶 POD 降 阶 误差 可 
小 于 0.9% ,可 代替 20 阶 的 未 降 阶 系统 。 

2) 在 获得 POD 降 阶 模型 时 ,可 基于 随机 激励 、 
扫 频 激励 和 谐 波 激励 啊 应 获取 ,而 脉冲 激励 所 得 到 
得 POD 降 阶 模型 无 法 得 到 符合 误差 标准 的 精度 。 
特别 在 使 用 谐 波 激励 响应 时 ,为 得 到 某 频率 的 最 佳 
降 阶 模型 ,应 选择 其 对 应 谐 波 激励 频率 进行 计算 。 

3) 低频 采样 频率 和 短 时 采样 时 长 也 可 达到 较 高 
的 降 阶 精度 ,不 必 追 求 高 频 采 样 频率 和 超 长 的 采样 
时 长 。 对 于 本 研究 对 象 , 取 100 Hz 最 低 分 辨 采样 频 
率 和 50s 的 采样 时 长 便 可 达到 与 1 000 Hz 采样 频率 
和 1 000 s 采样 时 长 相同 的 效果 。 

4) 本 研究 给 出 对 含有 高 斯 白 噪声 的 响应 信号 进 
行 降 阶 的 过 程 , 对 于 本 研究 对 象 ,通过 对 信和 号 进行 降 
噪 和 提高 阶 数 两 步 便 可 得 到 最 高 仅 为 2% 的 降 阶 预 
测 误差 。 为 实际 工程 中 的 应 用 提供 了 一 种 可 行 方 
法 ,并 能 够 使 得 该 方法 更 好 地 应 用 于 实际 中 。 

综合 可 知 ,可 使 用 基于 POD 的 模型 降 阶 方法 ， 
得 到 的 含 局 部 非 线性 悬臂 梁 低 阶 模型 ,对 进行 非 线 
性 结构 的 振动 控制 咒 设 计 有 实际 意义 。 
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